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ABSTRACT 

When measuring the mass profile of any given cosmological structure through inter- 
nal kinematics, the distant background density is always ignored. This trick is often 
rcfcrccd to as the "Jeans Swindle". Without this trick a divergent term from the back- 
ground density renders the mass profile undefined, however, this trick has no formal 
justification. We show that when one includes the expansion of the Universe in the 
Jeans equation, a term appears which exactly cancels the divergent term from the 
background. We thereby establish a formal justification for using the Jeans Swindle. 
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1 INTRODUCTION 

A small overdensity in an otherwise infinite homogeneous 
gravitating system (like any cosmological structure in the 
Universe) is affected by a basic inconsistency, namely that 
such a system cannot be in equilibrium, and at the same 
time obey the Poisson's equation which relates the gravita- 
tional potential to the density distribution. A constant grav- 
itational potential leads, via Poisson's equation, to a zero 
density (Jeans 1929; Zeldovich & Novikov 1971). The usual 
way to overcome this inconsistency is to assume that the in- 
finite homogeneous system does not contribute to the grav- 
itational potential, meaning that the gravitational potential 
is sourced only by fluctuations to this uniform background 
density. This assumption is called Jeans Swindle (Binney 
& Tremaine 1987, 2008; Kiessling 1999, 2003; Ershkovich 
2011). Following Binney & Tremaine (1987) "it is a swindle 
because in general there is no formal justification for dis- 
carding the unperturbed gravitational field". It is vindicated 
by the right results it provides, but it is generally considered 
a limitation to the formalism. 

The Jeans Swindle has several applications. Here we fo- 
cus on the Jeans Swindle in the context of the Jeans analysis 
of internal kinematics, which for instance is relevant for stel- 
lar motions in dwarf galaxies and galaxy motions in galaxy 
clusters. The aim of this work is to explain the "swindle" 
through a clean derivation of the Jeans equation, including 
the crucial expansion of the Universe. 

The Jeans equation describes systems in equilibrium, 
and it is therefore used to model for example dark mat- 
ter halos inside the virial region, where they can be treated 
as equilibrated systems. Dark matter (DM) halos can be 
seen as a matter excesses over the mean matter density of 



the Universe. This constant background density is the main 
contribution to the density distribution at large distances 
from the halo center (Tavio et al. 2008). We show that the 
contributions from the background density, the cosmologi- 
cal constant and the Hubble expansion, cancel each other. 
When omitting the constant background density (the nor- 
mal "swindle" ) one is actually excluding it together with the 
contribution from the expansion of the Universe. Thus, once 
we take into account the expansion of the Universe and the 
presence of the cosmological constant, we no longer need to 
invoke the Jeans Swindle. 



2 JEANS SWINDLE IN THE JEANS 
EQUATION 

The dynamics of DM halos, modelled as spherical and sta- 
tionary systems of collisionless particles in equilibrium, is 
controlled by the Jeans equation (Binney 1980) 
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where <j r is the radial velocity dispersion, /3 = 1 — o-J/of the 
velocity anisotropy, p the density distribution of particles 
and $ the total gravitational potential. 

The potential gradient is given by Poisson's equation 
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In the simple case of an isotropic velocity distribution (/3 = 
0), the solution to the standard Jeans equation (1) for the 
radial velocity dispersion is (from Binney 1980) 
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Thus, the only quantity required for the calculation of the 
radial dispersion is the density distribution of the halo. DM- 
only cosmological N-body simulations indicate that a double 
slope profile provides a reasonable fit to the density pro- 
files of halos within the virial radius (Navarro et al. 1996; 
Kravtsov et al. 1998). Since the integration in equation (3) 
extends all the way to infinity, we need a correct description 
of the DM distribution beyond the virial radius. The correct 
asymptotic value should be the mean matter density of the 
Universe pb g , given by 



Pb g = Quipc = 



3 £1 m H 
SttG 



(4) 



where Qm the matter density parameter, p c is the critical 
density of the Universe and H — a/a is the Hubble constant 
(a being the scale factor of the Universe). 

Therefore, the double slope profile, reaching zero den- 
sity at large distances from the cluster center, does not give 
the right results for the external regions (Tavio et al. 2008) . 
As a first approximation, we can write the density as given 
by the sum of a term ph that reproduces the inner part of 
the halo distribution and the constant background density 
that affects the profile only at large radii 



p(r) = p h (r) + p bg . 



(5) 



As an example, we consider a finite mass density profile for 
a cluster-size halo, the Hernquist (1990) profile 1 



phM = 



Po 



r/r v (1 + r/r v ) 3 



(6) 



where r v is the virial radius and po is the characteristic den- 
sity, that can be written in terms of the virial overdensity 
A as 



po 
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(7) 



In Figure 1 we show the profile given by eq. (6) (black solid 
line) and the sum ph(r)+pb s (green dash-dot line), where the 
asymptotic value is pb g (red dashed line). In the calculation, 
we set Q M = 0.24, A = 100, H = 73kms _1 Mpc" 1 and we 
fix a by imposing that the mass given by the density in 
eq. (6) within the sphere of radius r v corresponds to the 
virial mass 
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In Figure 1 the density is in units of Ap c and the radius is 
in units of the virial radius. 

Once we insert the expression (5) in the integral (3) we 

get 




Figure 1. Black solid line: Hernquist profile given by equa- 
tion (6). Red dashed line: background density given by equa- 
tion (4). Green dash-dot line: the sum of the Hernquist profile 
and the constant background density. 

which diverges, since the background mass diverges at large 
radii. In Fig. 2 we plot the solution (9) for different up- 
per limits i? max in the integral: i? max = 200 r v (blue short- 
dashed line), -R m ax = 500 r v (green dash-dot line), i? ma x = 
1000 r v (magenta dash-dot-dot line). 

This clearly indicates that when integrating to infinity, 
the integral will diverge. 

The usual trick to avoid the divergence is to omit the 
contribution of the background density to the potential gra- 
dient (2), i.e. to set V<E>b g = (Ershkovich 2011; Binney & 
Tremaine 1987, 2008). Physically, this amounts to assume 
that the gravitational potential is sourced only by fluctua- 
tions to the uniform background density. For this require- 
ment to be consistent with the Poisson's equation (2), the 
constant pb g in equation (5) has to vanish. This assumption 
is called the Jeans Swindle. It has no justification other than 
to overcome a mathematical difficulty. 

The Jeans analysis has always been performed by dis- 
carding the unperturbed density, so that equation (3) is just 
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The black solid line in fig. 2 corresponds to the computation 
of equation (10), which is not divergent. 

With this approach, the Jeans equation reproduces the 
radial dispersion of DM halos from cosmological simulations 
in the region of equilibrium (Sanchis et al. 2004; Wojtak 
ct al. 2005; Cuesta et al. 2008), and gives finite masses for 
dwarf galaxies (Strigari et al. 2010; Salucci et al. 2012). 
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1 Here we discuss the simple Hernquist profile for academic rea- 
sons. Using any other finite mass structure would lead to the same 
conclusions. 



3 WHY THE JEANS SWINDLE WORKS 

We wish to replace the Jeans Swindle by a formally correct 
analysis. Thus, we keep the background density pb g and its 
contribution in the gravitational potential. The gradient of 
the potential due to pb g can be put in the following form: 



d<E>bg 4 fi m 2 

— = r G Phs r=— Hr, 



(11) 
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Figure 2. Radial velocity dispersion profiles. The black solid 
line corresponds to the standard Jeans solution (10). Solutions 
of equation (9) corresponding to different upper limits J? ma x in 
the integral arc shown: R ma x = 200 r v (blue short-dashed line), 
flmai = 500 r v (green dash-dot line), -Rmax = 1000 r v (ma- 
genta dash-dot-dot line). The red long-dashed line is the solution 
when including cosmological expansion and background density 
(eq. [20]). The radial velocity dispersion is in units of the virial 
velocity and the radius is in units of the virial radius. 

where we have used equation (4). However, in order to be 
consistent, we need to take into account all effects due to the 
underlying cosmology. In a universe with a non-zero cosmo- 
logical constant A, particles also feel a repulsive potential of 
the form (e.g. Peirani & de Freitas Pacheco 2006; Nandra 
et al. 2012) 

d$ A 
dr 

where we have used the relation 

1U = ^. (13) 
Introducing the deceleration parameter 

~Y 

we can rewrite the total contribution of the cosmology to 
the gravitational potential gradient as 

d$ bg d$ A 
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Moreover, the Universe is not static, but it is subject to 
the Hubble expansion. The equation (1) holds for structures 
which have achieved dynamical equilibrium. This means 
that the radial, longitudinal and azimuthal bulk motions are 
not taken into account in its derivation. When excluding all 
these bulk velocity terms, also the Hubble flow, which DM 
particles are subject to, is discarded. The Hubble velocity, 
vh — Hr, might be neglected in the very inner region, but 
for large radii, it becomes important. Since the integration 
in equation (3) is all way to infinity, the inclusion of vh will 
affect the result. 

When we include the terms involving the mean radial 
velocity, the Jeans equation becomes the more general for- 
mula (Falco et al. 2012) 

d(pa 2 ) 
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In the most general case, v r would be the sum of the Hubble 
velocity and a peculiar infall velocity. If we keep neglecting 
the mean radial peculiar motion of particles just outside the 
equilibrated cluster, the radial velocity corresponds to vh 
only 



v r (r,t) = H(t)r. 



(17) 



It is straightforward to calculate the additional term involv- 
ing the Hubble velocity 



_ dv r dvr 
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where we used 
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We can now write equation (3) including all these cosmolog- 
ical terms 
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We thus see that the term — q H 2 s\h given by the Hub- 
ble velocity (18) cancels exactly the term qH 2 s|b g ,A given 
by the potentials of the background density and the cosmo- 
logical constant. In this way, we recover the same result as 
applying the Jeans Swindle, and in the Jeans solution the 
total mass is again M{s) = Mh(s). 

Formally, there is still a minor difference between the 
two approaches: the density involved in eq. (20) is still given 
by Ph(s) + Pbg, where pbg is not zero but instead given by 
eq. (4). However, this time it does not lead to any divergence, 
because Mh(s)/s 2 goes rapidly to zero at large distances. 
This can be seen in figure (2), where the solution (20) is the 
red long-dashed line and it is the same result we obtain from 
equation (10) (black solid line). For larger radii, the addition 
of the background density in the density profile can affect 
the result slightly. However, in the outer regions, where the 
halos are no longer equilibrated, the standard Jeans equation 
is anyway not used to reproduce the radial velocity disper- 
sion. Instead the generalized Jeans equation in eq. (16) must 
be used, including the infall motion of galaxies. One could 
also improve on this minor difference, by not including pbg 
at all radii, but instead a different form which takes into 
account that near the halo, what used to be background is 
now instead a part of the halo (see Falco et al. (2012) for a 
thorough discussion of this last point). 



4 COMOVING COORDINATES 

The equations describing the particle distribution and mo- 
tion can be written in comoving coordinates (Peebles 1980). 
The physical coordinates r and comoving coordinates x are 
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related by the universal time-dependent expansion parame- 
ter a(t) 

r = a(t)x. (21) 

When changing variables from the physical space to the co- 
moving one, the Poisson's equation becomes (Peebles 1980) 

V 2 <t> = 4nGa 2 [p(x)- p hg ], (22) 

where the gradient is with respect to x, and pb g is the mean 
mass density and 4>(x) is the potential contributed by the 
overdensity p — pb g . Therefore, in this coordinate system, the 
particle motion is already described in terms of the depar- 
ture from the constant background, and the swindle is not 
required. As we expect, taking into account the cosmological 
expansion in the physical space leads to the same result as 
moving to the expanding space. Equation (1) would be cor- 
rect if we replaced r with x and 3> with <j>, and using p given 
by (5), namely it is the correct Jeans equation in comoving 
coordinates. 

We conclude that our work is consistent with the comov- 
ing frame analysis by Peebles (1980). This confirms that the 
Jeans Swindle corresponds to accounting for the expansion 
of the Universe. 

5 CONCLUSIONS 

We have demonstrated that the Jeans Swindle is not an ad 
hoc trick, but it is the result of correctly combining the mean 
matter density and the expansion of the Universe. The di- 
vergent term from the background density, which in a static 
universe would lead to a divergent dispersion profile, is ex- 
actly cancelled by a term from the expanding universe. We 
have shown that the dispersion profile measured when as- 
suming no background and a static universe, is the same as 
the dispersion profile when including both the background 
density and the expansion. This means that we have estab- 
lish a formal justification for using the Jeans Swindle. This 
result holds for radii smaller than roughly the virial radius. 
For larger radii one has to include the effect of infalling mat- 
ter, which is done through a generalized Jeans equation, as 
will be presented in a forthcoming article (Falco et al. 2012). 
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